A new mechanism of neutron star radiation 
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We find a new mechanism of neutron star radiation 
wherein radiation is produced by the stellar interior. The 
source of radiation is oscillating neutron vortices in the su- 
perfluid core of a rotating neutron star. Thermally excited 
helical waves of vortices generate fast magnetosonic waves at 
the stellar crust. Near the crust bottom such waves reduce 
to a collisionless zero sound in an electron liquid, while near 
the stellar surface they behave as electromagnetic waves in a 
medium. The magnetosonic waves propagate across the crust 
and transform into electromagnetic radiation at the surface. 
The vortex contribution has the spectral index a w —0.45 and 
can explain nonthermal radiation of middle-aged pulsars ob- 
served in infrared, optical and hard X-ray bands. Detection 
of vortex radiation allows direct determination of the core 
temperature. Comparing the theory with available spectra 
observations we find that the core temperature of the Vela 
pulsar is T ~ 8 x 10 8 K, while the core temperature of PSR 
B0656+14 exceeds 2 x 10 8 K. This is the first measurement 
of the temperature of a neutron star core. The temperature 
estimate rules out equation of states incorporating Bose con- 
densations of pions or kaons and quark matter in these ob- 
jects. In principle, zero sound can also be emitted by other 
mechanisms, rather than vortices, which opens a perspective 
of direct spectroscopic study of superdense matter in the neu- 
tron star interiors. 



Properties of matter at densities much larger than nu- 
clear are poorly known and constitute a challenging prob- 
lem of modern science. They have broad implications of 
great importance for cosmology, the early universe, its 
evolution, for compact stars and for laboratory physics of 
high-energy nuclear collisions. Present searches of mat- 
ter properties at high density take place in several arenas. 
One of them is investigation of neutron stars. 

Neutron stars (NSs) are compact objects of radius 
about 10km and the mass of the order of solar mass. 
Cores of NSs consist of superfluid neutrons with some 
(several per cent by particle number) admixture of pro- 
tons and electrons. More dense NSs also possess an inner 
core. Its radius may reach several kilometers and central 
density can be much larger then the nuclei density. Sev- 
eral hypotheses about composition and equation of state 
of the inner core are discussed in the literature. One of 
them is appearance of S- and A-hyperons. The second 
hypothesis assumes the appearance of pion or kaon con- 
densation. The third hypothesis predicts a phase transi- 
tion to strange quark matter composed of almost free u, 
d and s quarks with a small admixture of electrons. 



The cooling rate and temperature of the central part 
of a NS substantially depends on the properties of dense 
matter. In this paper we find that NS radiation can be 
produced in its interior, rather then at the surface. De- 
tection of such radiation in the A— ray band allows direct 
determination of the core temperature. We find that the 
core temperature of the Vela pulsar is T » 8 x 10 8 K, 
while the core temperature of PSR B0656+14 exceeds 
2 x 10 8 K. This estimate excludes exotic equation of states 
incorporating Bose condensations of pions or kaons and 
quark matter in these objects. A systematic study of 
pulsar radiation in the A— ray band can locate objects 
with fast cooling core which could be candidates for NSs 
with exotic states of matter. 

Zero sound in electron liquid — For middle aged NSs 
the effective frequencies of electron-electron and electron- 
phonon collisions in the stellar crust are less than 10 12 Hz. 
In the core the collisions of electrons with baryons are 
suppressed by neutron superfluidity and proton super- 
conductivity. As a result, at large frequencies the elec- 
trons in the stellar interior can be treated as an indepen- 
dent system of Fermi particles which move in collision- 
less regime. In such regime a zero sound is a possible 
collective excitation of Fermi liquid Speed of zero 
sound u s is larger than the electron Fermi velocity Vf- 
In the presence of magnetic field H equations of elec- 
tron motion reduce to equations of magnetic hydrody- 
namics in which u s enters the equations instead of speed 
of usual sound and the Alfven velocity ua is estimated 
in terms of the electron density p e : ua = H / ' \J \np e . 
Magnetic field at the stellar crust couples electron sound 
and electromagnetic wave into one excitation known as 
a fast magnetosonic wave (FMW). Near the crust bot- 
tom ua <C u s rs c, where c is the speed of light, and the 
FMWs reduce to zero sound. However, near the stellar 
surface (at small density) ua > u s and the FMWs behave 
as electromagnetic waves in a medium with dielectric con- 
stant e = l + c 2 /u A . As a result, at the stellar surface the 
FMWs transform into electromagnetic radiation in vac- 
uum the same way as refraction of usual electromagnetic 
waves at the dielectric-vacuum interface. Typically near 
the stellar surface ua 3> c and, hence, the transmission 
coefficient is close to 1. 

The key idea of the paper is that NS interior is trans- 
parent for zero sound (FMWs). We assume that the crust 
matter is in a crystal phase: T < T m = 2.5 x 10 7 pg^ 3 K, 
where T m is the melting temperature, p§ is the matter 



density in units 10 g/cm 



Near the stellar surface 
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there is a phase transition between a gaseous atmosphere 
and a condensed metallic phase and the atmosphere has 
negligible optical depth H|. The density of the metal- 
lic phase at the surface depends on magnetic field. E.g., 
for iron crust and H = 10 12 Gs the surface density is 
p « 10 3 g/cm 3 §. 

For a weakly interacting Fermi gas, such as a gas of 
electrons in the NS interior, the attenuation of zero sound 
is exponentially small ||. However, in the presence of 
magnetic field the FMWs are accompanied by oscillation 
of electric field which results in partial damping of FMWs 
near the stellar surface. The attenuation increases with 
decreasing the matter density. However, at p§ < QSHf^ 
the system as a whole occupies only a small number of 
Landau levels 0. At such densities the electron mo- 
tion becomes effectively one dimensional and the elec- 
tron scattering, as well as magnetosonic attenuation, is 
suppressed by magnetic field. The maximum damping 
occurs in a boundary region where magnetic field still 
does not suppress the electron scattering. FMWs pass 
the critical region and reach the surface without attenu- 
ation if in the region 

|sin0| < 0.004iJ 2 2 . 9l f /T 6 5/2 , (1) 

here 6 is the angle between the wavevector and the mag- 
netic field, <7i5 is the acceleration of gravity in units 
10 15 cm/s 2 f§. For T = 10 6 K, H = 4 x 10 12 Gs and 
g = 10 15 cm/s 2 only waves with \9\ < 3.7° reach the 
stellar surface. Below we show that zero sound in NS 
interior is generated at the core-crust interface by super- 
fluid neutron vortices which undergo thermal helical os- 
cillations. Vortices themselves are excited by scattering 
of single-particle excitations on the vortex cores. Sound 
waves propagate across the stellar crust towards the sur- 
face. Due to anisotropic attenuation in the surface region 
and star rotation the observed vortex radiation is pulsed. 
According to Eq. (Q) , for colder NSs the shape of vortex 
pulses becomes broader. 

Vortices in neutron stars — If a NS rotates with an 
angular velocity Q =Qz a periodic vortex lattice forms 
in a superfluid neutron phase of the stellar core. The 
characteristic size of the vortex is £ = 10 -12 cm and 
the total number of vortices for fl/2ir = 10s _1 is N = 
2m n R 2 Vt/h w 10 17 , where m n is the neutron mass and 
R w 10km is the radius of the superfluid. Let us con- 
sider a single straight neutron vortex line with length L. 
Normal modes of the vortex are helical waves Q 

r = Aexp(iuit + ik\\z), oj = Tiki ln(l/fc||£)/4m„, (2) 

where A is the amplitude and k\\ is the wave number. 
Eqs. (0) are valid for fc||£ -C 1, or lu <C lu c — n/4m n £ 2 
10 20 Hz and describe rotation of a helix with the angular 
frequency lu. At distances r < £ (we use cylindrical coor- 
dinates r, <fr, z) the helical vortex motion can be described 
in terms of the perturbation in the velocity potential $' 



$' = i Ah exp [i(u>t + k\\z — </>)] /2m n r. (3) 

In the bulk the function <E>' should satisfy a linearized 
equation of superfluid hydrodynamics, which is, in fact, 
the wave equation of sound propagation 

d 2 $'/dt 2 - c 2 V 2 $' = 0, (4) 

where c s is the speed of sound in superfluid. In the limit 
lu -C lu c one can omit the first term in Eq. ([|). Then so- 
lution of this equation which satisfies the boundary con- 
dition ([}]) has the form 

= iAhk\\Ki{k\\r) exp [i(wt + k\\z - 4>)] /2m n , (5) 

where K\ is the modified Bessel function. The solution 
exponentially decreases at r 3> l/fc||, that is perturbation 
in the superfluid velocity V = V$ is localized near the 
vortex core and no energy is radiated. 

Generation of zero sound at the core boundary — The 
bulk of superconducting protons in NS core do not ro- 
tate by forming vortices |8|. However, due to interac- 
tions between protons and neutrons, neutron superfluid 
velocity generates a superfluid current of protons in the 
vicinity of neutron vortices (drag effect). At the interface 
between the stellar crust and the outer core the matter 
undergoes a first order transition from an Aen phase to 
a uniform liquid of neutrons, protons and electrons (npe 
phase). Nuclei in the Aen phase form a Coulomb crystal. 
In the npe phase the proton liquid is superconducting, 
while in the Aen phase protons constitute a part of nu- 
clei. According to Eq. (||), helical vortex motion does 
not generate sound waves inside the stellar core. How- 
ever, such waves are generated at the interface between 
the superconducting and Aen phases. 

Due to the drag effect helical waves of a neutron vor- 
tex produce density oscillation of protons in the super- 
conducting phase. Free electrons screen the electric field 
and oscillate together with protons. In Aen phase there 
is no drag effect and electron motion is not coupled with 
neutron vortices. To estimate the power of radiated zero 
sound we note that kinetic energy of ultrarelativistic elec- 
trons at the crust bottom is much larger than their in- 
teraction energy. Hence, the radiated power is approxi- 
mately the same as in the limiting case of non interacting 
electrons. In this limit, the electrons move in a ballistic 
regime and the energy flux across the crust-core interface 
is given by the surface integral Q = J E e vpdS, where 
vf ~ c and E e is the contribution to the electron energy 
density due to the helical vortex motion. Density oscilla- 
tion of electrons (which is accompanied by large change 
in their Fermi energy) determine the helix energy in the 
stellar core which is equal to J E e LdS. As a result, the 
average energy flux of sound waves produced by a helical 
wave of a single vortex is 

Q = lJ EeLdS = I Mto/faq - 1] ' (6) 
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where T is the core temperature. One should note that 
the wave length of sound is much larger than the size 
of the area from which it is generated and, hence, the 
sound waves have approximately spherical front. Taking 
into account that the number of vortex modes within 
the interval dku is Ldk\\/2n and the number of vortices 
N = 2rn n R 2 n/ti, we obtain the spectral density of sound 
waves power radiated from the oscillating vortex lattice 
in a half space 



PvM = 



x/2, 



cm i. 



nh 3/2 y/hln(u) c /u;) [exp(huj / k B T) - 1] 
Fig. 1 shows vortex mechanism of NS radiation. 
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FIG. 1. Mechanism of a NS radiation. Thermaly excited 
helical waves of neutron vortices in the superfluid core pro- 
duce FMWs in the stellar crust. FMWs propagate across 
the crust and transform into electromagnetic radiation at the 
star surface. Mainly the radiation comes out from regions 
with strong magnetic field (near magnetic poles). 

One should note that electromagnetic radiation pro- 
duced by vortices can not exceed black body radiation 
produced by the surface 4irR 2 with the same temperature 
T (Kirchhoff's law). We assume that vortices are excited 
thermally and if at some frequency the vortex radiation 
becomes comparable with radiation of a black body this 
means that at such frequency the rate of thermal excita- 
tion imposes the main restriction on the radiated power: 
vortices radiate maximum power which can be pumped 
from the thermal reservoir. Kirchhoff's law results in the 
following limitation at which Eq. (Q) describes vortex 
radiation 



w/2vr > 0.24h- 3/5 c 6/5 m 3 J 5 n 2 / 5 \n- 1/5 {u c /uj). (8) 

For Q./2-k = 4s" 1 we obtain uj/2it > 1.5 x 10 14 Hz. At 
lower frequencies the spectrum of vortex radiation follows 
Planck's formula. 

To compare our theory with observations it is conve- 
nient to represent the spectral density of vortex radiation 
P v (lo) as a power law P v (u>) cx uj a , where a is the spec- 
tral index. The y/ln(u) c /uj) function in Eq. (^) shifts 
the spectral index by a small value 1/2 hx(u> c /cJ). The 
shift depends weakly on u> and changes the spectral in- 
dex from a ~ —0.46 in the optical band to a ~ —0.43 in 
the A— ray band. Apart from vortex contribution there 
is thermal radiation from the NS surface. We estimate 
thermal radiation assuming helium atmosphere and use 
numerical results obtained by Romani |J. The model 
fits well the thermal spectrum with only one parameter, 
the effective surface temperature T c g. The total radia- 
tion from a NS in unit solid angle is given by the sum of 
vortex and thermal Pth components: 



P(«) =aP v (w)/27r + Pth(w). 



(9) 



Here we introduced a dimensionless free parameter a rs 
0.1 -j- 1 which takes into account partial absorption of the 
vortex radiation near the stellar surface and geometri- 
cal effects related to unknown magnetic field distribution 
and position of the line of site. 

Discussion — Fig. 2 compares the observed radia- 
tion spectrum of middle-aged pulsars PSR B0656+14 
(n/2ir = 2.6s- 1 ) and Vela {fl/2ir = 11.2s- 1 ) with the 
spectrum predicted by our theory. Thermal radiation of 
the stellar surface dominates in the ultraviolet and soft 
A— ray bands. The effective surface temperature is T e g = 
4.9 x 10 5 K for PSR B0656+14 and T cS = 7.8 x 10 5 K for 
Vela. The NS radius R s (in km) is related to the pul- 
sar distance D (in kpc) as R s = 53\/l + zD for PSR 
B0656+14 and R s = 41y/T+lD for Vela (z is the NS 
redshift). The vortex contribution (dash line) dominates 
in infrared, optical and hard A— ray bands. In the far 
infrared band the radiation spectrum changes its behav- 
ior and follows Planck's formula. The sum of the vor- 
tex and thermal components is displayed with the solid 
line. For PSR B0656+14 we take the core temperature 
T = 6.4 x 10 8 K and the geometrical factor a = 0.18, while 
for Vela T = 8 x 10 8 K and a = 0.33. The observed broad- 
band spectrum is consistent with our model for typical 
NS parameters, which suggests that the vortex mecha- 
nism of radiation operates in a broad frequency range 
from IR to hard A— rays. Also in the optical band the 
radiation of middle-aged pulsars was found to be highly 
pulsed with pulses similar to those in the hard A— ray 
band JTo| , |TT[ | . This agrees with the vortex mechanism 
which predicts pulsed radiation. 
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FIG. 2. Broadband spectrum of PSR B0656+14 [12] and 
the Vela pulsar [13,14]. Solid line is the fit by the sum of the 
vortex and thermal components. Thermal radiation of the 
stellar surface dominates in the ultraviolet and soft X— ray 
bands, while the vortex contribution (dash line) prevails in 
infrared, optical and hard X— ray bands, where its spectrum 
has a slope a ~ —0.45. In the far infrared band the spec- 
trum changes its behavior and follows Planck's formula with 
P tx lo 2 . 

It is worth to note that the spectral index of vortex 
radiation is a fixed parameter in our theory. If a single 
mechanism of nonthermal radiation operates in the broad 
range, then the continuation of the X— ray fit to the op- 
tical range determines the spectral index of the nonther- 
mal component with a big accuracy. Good quantitative 
agreement of our theory with the observed spectral in- 
dex serves as a strong evidence that the vortex mecha- 
nism is responsible for radiation of middle-aged NSs in 
IR, optical and hard X— ray bands. The vortex contribu- 
tion exponentially decreases at fiui > ksT. Observation 
of such spectrum behavior allows direct determination of 
the core temperature T. For PSR B0656+14 and the Vela 
pulsars the power law spectrum in the hard X— ray band 
shows no changes up to the highest frequency 2 x 10 18 Hz 
at which the data are available (see Fig. 2). This in- 
dicates that temperature of the NSs core is larger than 
2 x 10 8 K. Measurements of spectra of middle-aged pulsars 
in the 10 18 — 10 19 Hz range are needed to search for pos- 
sible manifestations of the core temperature. Recent ob- 
servation of the Vela pulsar with RXTE has covered the 
frequency band (0.49 — 7.3) x 10 18 Hz. These data in com- 
bination with OSSE observations, (1.7-13.8) x 10 19 Hz, 
allows us to estimate the temperature of the Vela core. 
Light curves of the Vela pulsar have several peaks. We as- 
sociate the vortex radiation with the second optical peak 
(see Fig. 1 in Q and Fig. 4 in Q). To estimate the 
core temperature one should trace the evolution of the 
vortex peak in the energy interval IkeV— lMeV. If we as- 



sume that the amplitude of the vortex peak follows Eq. 
([?]), then we can find total (integrated over frequency) 
radiation intensities of vortices in different energy bands. 
Then we compare the detected integrated intensities of 
the vortex peak in the Rossi and OSSE bands |ll| with 
those predicted by Eq. (R). We found that if there was 
no temperature decay of the vortex peak the total vor- 
tex intensity in the OSSE band should be 3 times larger 
than those actually observed. This indicates on the de- 
cay of the vortex spectrum with the core temperature 
r«8x 10 8 K. 

The estimate of the core temperature allows us to make 
a conclusion about the interior constitution of the NSs. A 
hot NS cools mainly via neutrino emission from its core. 
Neutrino emission rates, and hence the core temperature, 
depend on the properties of dense matter. If the direct 
Urea process for nucleons is allowed such NS cools to 
10 8 K in weeks J!?]]. However, the characteristic age of 
PSR B0656+14 is about 10 5 yre, while the Vela pulsar 
is about 10 4 years old. So, the NS core cools down to 
10 8 K at least 10 6 times slower then the rate predicted by 
the direct Urea process. This estimate excludes equation 
of states incorporating Bose condensations of pions or 
kaons and quark matter. The point is that the neutrino 
emission processes for these states may be regarded as 
variants of the direct Urea process for nucleons Jl7| . As 
a result, all these states give rise to neutrino emission, 
though generally smaller, comparable to that from the 
direct Urea process for nucleons which inconsistent with 
the discrepancy of the cooling rate in the factor 10 6 . 

Detection of vortex radiation opens a possibility to 
study composition of NS crust. Since FMWs gener- 
ated by vortices propagate through the stellar interior 
the spectrum of vortex radiation should contain (red- 
shifted) absorption lines which correspond to low en- 
ergy excitation of nuclei that form NS crust. E. g., 
the 57 Fe nucleus has an excited state with the energy 
14.4kev= 3.5 x 10 18 Hz, which would produce an absorp- 
tion line in the hard A— ray band if the core temper- 
ature is greater than 1.7 x 10 8 K. Bottom layers of the 
crust may contain exotic nuclei with the mass number 
up to 600 and the core radiation creates a perspective 
to study their properties. Another challenging problem 
is interaction of zero sound with exotic states of matter. 
The point is that generation of zero sound by vortices is 
only one of the possible mechanisms. If an exotic state 
of matter absorbs zero sound at some frequency, then, 
according to Kirchhoff's law, it will radiate zero sound 
at this frequency. As a result, spectrum of stellar ra- 
diation must contain characteristic emission lines corre- 
sponding to such processes which opens a perspective of 
direct spectroscopic study of superdense matter. A de- 
tailed version of our theory is available in |(| . 
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